function L = BB98d_loss_time_domain(x,T_a,dT_a_dt,T_o,l_use_loss,P)

    PP = BB98d_predict(T_a,dT_a_dt,x);

    if P.do_season == 1
        intv = 12;
    else
        intv = 1;
    end
    PP.T_o_m = CDC_demean(PP.T_o_m,1,intv);
    T_o      = CDC_demean(T_o,1,intv);
    
    
    dif      = PP.T_o_m - T_o;
    
    l_use_loss(1:36)     = 0;
    dif(l_use_loss == 0) = nan;
    dif                  = dif - nanmean(dif);
    L                    = nanmean(dif.^2) * 100;
end
